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Abstract 

We study Andreev bound states and vr-junction transition in a supercon- 
ductor / quantum-dot / superconductor (S-QD-S) system by Green function 
method. We derive an equation to describe the Andreev bound states in 
S-QD-S system, and provide a unified understanding of the 7r-junction transi- 
tion caused by three different mechanisms: (1) Zeenian splitting. For QD with 
two spin levels and , we find that the surface of the Josephson current 
I{(f> = ^) vs the configuration of {E^,Ei) exhibits interesting profile: a sharp 
peak around E^ = E^ = 0; a positive ridge in the region of E-^ ■ E^ > 0; and a 
negative, flat, shallow plain in the region of Ey E^ < 0. (2) Intra-dot interac- 
tion. We deal with the intra-dot Coulomb interaction by Hartree-Fock approx- 
imation, and find that the system behaves as a vr-junction when QD becomes a 
magnetic dot due to the interaction. The conditions for vr-junction transition 
are also discussed. (3) N on- equilibrium distribution. We replace the Fermi 
distribution f{u)) by a non-equilibrium one ^ [f{uj — Vc) + /(w + Vc)], and al- 
low Zeeman splitting in QD where E^ = —E^ = h. The curves of = ^) 
vs Vc show the novel effect of interplay of non-equilibrium distribution with 
magnetization in QD. 
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PACS numbers: 74.50.+r, 73.40.Gk, 73.20.Dx, 72.15.Nj. 



I. INTRODUCTION 

Superconductivity has the nature of quantum condensation in a macroscopic scale, which 
can be described by a wavefunction with phase factor e''^ . When two superconductors are 
weak hnked, the phase difference will manifest in the dc Josephson current, with the current- 
phase relation / = /csin(0i — 02). If the weak link area is controlled by certain external 
conditions, the magnitude of the critical current Ic may be either suppressed or en- 
hanced 0. In some occasions, even the sign of Ic may be reversed 0, or equivalently the 
phase factor sin(0i — 02) changes to sin(0i — 02 + ^r), referred to as the vr-junction tran- 
sition. One of the simplest example is the so-called superconducting quantum diffraction, 
in which a superconductor / insulator / superconductor tunnel junction is tuned by an 
external magnetic field. The dc Josephson current vs the magnetic flux $b has the form 
I = h sin(0i — 02), and changes its sign at every ^b/^o equal to an integer. 

Unfortunately, the vr-junction transition cannot be directly detected in the two terminal tun- 
nel junction because the current source is used in the measurement to control supercurrent 
rather than phase difference. However, in a mesoscopic superconductor / normal-metal / 
superconductor (SNS) junction with the N-region coupled to normal electrode(s) 0J§], the 
phase difference can be determined independently by the coherent Andreev reflection cur- 
rent which is proportional to je^"^^ + e^'^^j^ = 2 [1 + cos(0i — 02)]- Thus, direct observation 
of the vr-junction transition becomes accessible. 

In mesoscopic SNS junctions, supercurrent is conducted through the N-region by Andreev 
reflection (AR) process 0. The energy gaps of two superconducting electrodes serve as two 
"mirrors", reflecting electron into hole and hole into electron. For ballistic SNS junctions, 
discrete Andreev bound states are formed in the N-region, each state carries positive or 
negative supercurrent. For diffusive SNS junctions, the so-called current carrying density of 
states (CCDOS) plays the similar role, which also has positive and negative contributions 
to the supercurrent. Recent experiment demonstrated the vr-junction transition in the 
diffusive SNS junction by applying control voltage on the N-region. In fact, the biased 
normal reservoirs across mesoscopic N-region induce a non-equilibrium distribution in the 
N-region, and make the occupied fraction of CCDOS deviating from the equilibrium one. 
When the control voltage is excess a certain value, the non-equilibrium distribution has so 
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much weight on the negative part of CCDOS that the total current reverses its sign. Many 
theoretical works have been addressed on this issue, either for ballistic SNS junctions Ij^-Q, 
or for diffusive SNS junctions [11|-|T^. 

In addition to non-equilibrium distribution, there is a completely different mechanism to 
realize the vr-junction transition, i.e., coupling superconductors by an Anderson impurity or 
an interacting quantum dot (QD). The works of Glazman et al. [|I^] and Spivak et al. 
revealed that when the impurity is single occupied, the sign of Josephson current for infinite 
Coulomb repulsion is opposite to that without the repulsion. Ishizaka et al. obtained 
the condition for vr-junction transition, by using non-crossing approximation and varying 
the strength of the Coulomb repulsion, the bare level position, the tunneling strength, and 
the temperature. Rozhkov et al. [|18[ analyzed the system in a non-perturbative way, and 
found a novel intermediate phase in which one of = and = vr is stable while the other 



is metastable, with the energy E{(j)) having a kink somewhere in between. Clerk et al. |T9 



studied the case of infinite-?/ and the regime where the superconducting gap A and the 
Kondo temperature Tx are comparable, showing that the position of the sub-gap resonance 
in the impurity spectral function develops a strong anomalous phase dependence, and the 
TT-junction behavior is lost as the position of the bound state moves above the Fermi energy. 

Recently, there are growing interests in the physics of superconductor in contact with 
ferromagnetic material, and the following works revealed another approach to achieve vr- 
junction. Prokic et al. |^ presented a theory of the vr-junction transition in atomic-scale 
superconductor / ferromagnet (S/F) superlattices. They found that the critical Josephson 
current has a non-monotonic dependence on the exchange field h in the ferromagnetic layer, 
becoming zero at the critical value, corresponding to the transition between = and 
= vr in the ground state. Yip et al. [^] and Heikkila et al. |22] demonstrated that the 
supercurrent through a mesoscopic SFS junction oscillates with an exponential decreasing 
envelope as a function of the exchange field or the distance between the electrodes. They 
also proposed that the suppressed supercurrent by the exchange field can be recovered by a 
proper non-equilibrium distribution. 

With these in mind, we are curious whether there is any relationship among the above 
three mechanisms for vr-junction transition. Motivated by this question, we investigate the 
following cases of vr-junction transition in a superconductor / quantum dot / superconductor 
(S-QD-S) system and provide a unified picture based on Andreev bound states. In section 
II, we study the vr-junction transition caused by Zeeman splitting. Assuming the QD with 
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two spin levels and Ei, we find that the surface of the Josephson current /(f) vs the 
configuration of {E-^,Ei) exhibits interesting profile: a sharp peak around E-^ = E^ = 0; 
a positive ridge in the region of E-^ ■ E^ > 0; and a negative, flat, shallow plain in the 
region of E-^ ■ E^ < 0. In section III, we study the vr-junction transition caused by intra- 
dot interaction. We model QD by Hdot = Eq X^o-'^tco- + Uriirii, and handle the interaction 
term by Hartree-Fock approximation. Thus this case is reduced to the first one except a 
self-consistent calculation for {n„). We show that the vr-junction transition occurs when QD 
becomes a magnetic dot due to the interaction. The conditions for vr-junction transition are 
also discussed. In section IV, we study the vr-junction transition caused by non-equilibrium 
distribution in QD. By replacing the Fermi distribution f{uj) with a non-equilibrium one 
i [f{uj — Vc) + f{uj + V^)], and allowing Zeeman splitting in QD as E^ = —E^ = h, we find 
that for h = 0, the supercurrent reverses its sign when control voltage Vc is excess a certain 



value, which is agree with previous work [|2J] and the experiment 0. For h ^ 0, the curves 
of Josephson current vs control voltage show a novel effect of interplay of non-equilibrium 
distribution with the magnetization in QD. Finally, we summarize our understanding of 
vr-junction transition in S-QD-S system in section V. 



II. QD WITH TWO SPIN LEVELS 



A. model Hamiltonian and formulation 



In this section, we study the S-QD-S system modeled by the following Hamiltonian: 

H = HL + Hn + Hdot + HT , (1) 

ka k 

Hr = Y. 'v^lX^ + E [^^~"^%b^-r>i + h.c 

pa p 

Ht = Y \^La\„Ca + /i-c] + Y VR^laC<^ + h.c] , 

ka pa 

where Hi, Hji describe the left and right superconducting leads with phase difference (pi — 'pR 
2^, Hdot describes the quantum dot with two spin levels, and H^ is the coupling between 



the quantum dot and the superconducting leads. 
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Since Josephson current can be expressed in terms of the Green functions of the QD, 
we first derive the Green function by solving Dyson equation. Following the formulation 
in 



p^ , we denote G and g as the Green functions of QD in Nambu representation, with 
and without the coupling to leads, respectively; and denote S as the self-energy due to the 
coupling between QD and the leads. The retarded Green function (Fourier transformed) of 
isolated QD is 

\ 











(2) 

The retarded self-energy (Fourier transformed) under the wide bandwidth approximation 



can be derived as [E3 



Sl/^(^) 



\ 



(3) 



where F^,//? is the coupling strength between the superconducting leads and QD, defined by 
Vlir = 2iTNL/jit'^^j^, in which Nl and Nr are the density of states in the left and right leads 
in normal state. The factor p{uj) is defined as. 



p{uj) 



\uj\ > A 
\uj\ < A 



(4) 



Notice that p{uj) is the ordinary dimensionless BCS density of states when |ci;| > A, but 
has an imaginary part when |u;| < A, corresponding to Andreev reflection process within 
the superconducting gap. For simplicity, we assume that the two superconducting leads are 
identical except a phase difference . Let 0i = |, 0/? = — ^, Ff, = rR = r, then we obtain 
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By using Dyson equation, the retarded Green function of QD can be obtained as, 

\ 
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(7) 



where A = A{u) defined as 
A{uj) = det 



(.922 ^22)(5'll ^11, 



^12 ^21 



(8) 



The general current formula for a mesoscopic hybrid multi-terminal system has been 
derived in . For the time- independent case and QD with two spin levels under consider- 
ation, the current formula can be rewritten in a compact form (in units of e = h = 1): 



where [CD] _ ^ ^ , ^ ^ , ^ iu-22 
transformed 2x2 Nambu matrices. Since I 
current formula can be further reduced to 



/ ^2Re[GS,/^]<_^^ , (9) 
[ ]22, and G, Sl/r are the Fourier 



Jii 
h 



-Ir in the stationary transport, the 



I = -{Il-Ir) 
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(10) 



with S 



'R 



Applying the fluctuation-dissipation theorem, one has 



G<=f{u) [G^-G"] and S< = f{uj) 
distribution function. Notice that (G' 
can be simplified to 



where fiuj) = l/(e'^'^ + 1) is the Fermi 



G", (S^ 



S , the expression in the integrand 



Re 



GS 



f{uj)2 sin0- 



r2^2 



11-22 ' ' ' — A"^ 

Consequently, the Josephson current is expressed as 



— Im 



Aiuj) 



(11) 



/duj 



(12) 



in which the current carrying density of states (CCDOS) j{uj) is defined by 



r2^2 



- A2 



Im 



Aiu) 



(13) 



Because the singularities of j{uj) lie in the same half-plain, CCDOS j{uj) satisfies the condi- 
tion / i{uj)duj = 0. 

Since ^^{uj) is purely imaginary when > A while purely real when \uj\ < A, so A{uj) 
has finite imaginary part when \uj\ > A, while infinitesimal imaginary part when Ic^l < A. 
Correspondingly, the Josephson current can be divided into two parts, contributed from the 
continuous spectrum and from the discrete spectrum, respectively: 
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I = Ic + Id , 

Ic = '2 sin I 
Id = 2 sin I 



(14) 



27r' 



2^ 



We shall show in the appendix that when |ti;| < A and 7^ 0, the equation A{lj) = has 
two real roots denoted by -Ei = A sin 6'iand E2 = A sin 62, where 61 and 62 are the two roots 
of the following equation, 



sin 6' + — tan 61 - ^ j ( sin 6* + — tan 61 + ^ j cos^ ^ ~ ^ 2 ^ 



(15) 



with 6 G (— f , f )• Equation (15) completely determines the properties of Andreev bound 
states. The roots i^iand E2 are just the Andreev bound states. In the range of |ci;| < A, 



A{uj) can be written as a{uj){uj — Ei+iO^){uj — E2+i0~^), and 



Im -rrr 



is reduced to 



Im 



1 



A{u;) 



IT 
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^ 5{uj - ^1) + — ^-^5{uj - E2) 



A'{E2 



Finally, the Josephson current through S-QD-S system can be expressed as. 



I = Ic + h 

Ic = 2 sin 



(16) 



(17) 



smi 
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Im 
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2^ 



E2-A2A'(^i) El- A^ A{E2)_ 

This current formula will be used in the following numerical study [BO]. 



B. numerical results and discussions 

Now we discuss the numerical results for Andreev bound states and the Josephson cur- 
rent. In all numerical studies of this paper, we take e = /i = fc^ = 1, set A = 1, i.e., measure 
all energies in units of A, let F = 0.1 for symmetric and weak coupling case, and fix the 
phase difference = f . 

Fig.l presents the solution of Eq.(15), i.e., Andreev bound states of S-QD-S system. 
In the limit of F — > 00, Eq.(15) gives Ei = |cos|| and E2 = — |cos||, the well-known 
Andreev bound states for a clean superconducting point contact [|TU|. Conversely, in the 



limit of F ^ 0, the two roots of Eq.(15) are Ei = E^ and Ei = —E^, i.e., the bare levels 
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of electron with spinj and hole with spinj of the QD. (We describe spinf quasiparticle 
in electron language and spin| quasiparticle in hole language due to the choice of Nambu 
representation.) For the case of F <^ A under consideration, where QD is weakly coupled 
with the superconducting leads, the solutions of Eq.(15) depend strongly on the configuration 
of QD levels {E^, Ei), but weakly on the phase difference 0. The surfaces of Ei vs {E^, Ei) 
and E2 vs {E^, Ei) for </> = f are shown in Fig. la and Fig. lb, respectively. In this case, the 
electron level E-^ and the hole level —Ei are coupled by AR tunneling. Therefore, Andreev 
bound states can be viewed as hybrid states of E^ and — -Ej,, and an energy gap of the order 
F is opened where E-^ and —E^ are equal (see Fig.lc and Fig. Id). Further study on the 
relations of Ei vs (p and E2 vs (p can provide the information of supercurrent carried by each 
Andreev bound state (not shown). 

Fig. 2 presents the surface of Josephson current / vs the configuration of [E-^, E^). First, 
the surface is symmetric to the diagonal lines of E-^ = E^ and E-^ = —E^, which reflects the 
symmetry between electron and hole and the fact that supercurrent are non-spin-polarized. 
Second, the surface exhibits interesting profile: a sharp peak around E^ = E^ = 0; a. 
positive ridge in the region of E-^ ■ E)^ > 0; and a negative, fiat, shallow plain in the region 
of E-^ ■ E^ < 0. Note that the positive ridge and negative plain share a sharp edge. Third, 
consider two special but typical cases, E-^ = E^ = Eq (Fig.2b), and E^ = —E^ = Eq (Fig. 2c), 
which are actually the diagonal cuts of Fig. 2a. For the case of E-^ = E^, the I vs Eq curve 
has a peak at E'q = with a width of F, which reproduces the result for QD with one spin- 
degenerate level [|^. For the case of E-^ = —E^, the / vs Eq curve has the same maximum 
at Eq = 0, but jumps suddenly from the positive maximum to a small negative value around 
Eq = ±r, which is quite similar to the curve of the critical current vs the exchange field in 



an atomic scale S/F superlattice (see Fig.2 of |20|). 

To understand these features, we plot CCDOS j{u!) for different configurations of {E-^, E^) 
in Fig.3. The supercurrent can be expressed as / = 2sin0j ^f(uj)j(uj), and only the 
spectrum of u; < devotes to the current at zero temperature. For the spectrum of E^ = 
Ei = , jiuj) has two 5-function type discrete spectrum within the superconducting gap, 
corresponding to two Andreev bound states. They carry supercurrent with opposite signs, 
positive for Ei < 0, and negative for E2 > 0. j{uj) also has continuous spectrum outside 
the superconducting gap, negative and positive ior u < — A and u > A, respectively. Since 
the contribution from the discrete spectrum is much larger than that from the continuous 
one, the current peak at E^ = Ei = is mostly contributed from Ei. For the spectrum 



8 



of = El = Eq ^ 0, El and E2 move toward ±A, symmetric to the Fermi surface 
(see also Fig.lc). The contribution from the discrete spectrum Ei decreases continuously 
with Eq, corresponding to the F-width broadening of the Eq = Q peak in Fig.2c. For the 
spectrum of E^ = —E^ = Eq 0, however, Ei and E2 move in the same direction. When 
Eq < — F or Eq > F, both Ei and E2 are below or above the Fermi surface (see also Fig. Id). 
As a consequence, they have little net contribution to the supercurrent, and the relatively 
small negative continuous spectrum of < —A dominates. Because the crossover of Ei 
and E2 from different sides of the Fermi surface to one side occurs abruptly, sudden jumps 
between the positive maximum and the negative valleys appear in Fig.2d. Similarly, one can 
understand the whole surface of Fig. 2 with the help of the Andreev bound states in Fig.l 
and the properties of CCDOS in Fig.3. 

The above results are for zero temperature, the temperature dependence of the Josephson 
current is shown in Fig. 4. Notice that the supercurrent is sensitive to the temperature. The 
height of E^ = E^ = peak (referred as Iq) decreases rapidly with the increase of the 
temperature, Jq = 0.06, 0.004, 0.001 for T = 0, 0.25, 0.50, respectively. The sharp edge 
between / > and / < is also smeared out at finite temperatures. This reason is, at finite 
temperature not only the Andreev bound state below the Fermi surface but also the one 
above the Fermi surface has contribution to the supercurrent. 



III. QD WITH INTRA-DOT INTERACTION 

A. Hartree-Fock approximation 

In this section, we investigate the vr-junction transition caused by the intra-dot in- 
teraction. The quantum dot with Coulomb interaction can be described by H^ot = 
Eq^^cIc^ + Un^ni. As in the problem of local moment in nonmagnetic metals [^, we 



deal with the interaction term by Hartree-Fock approximation (HFA), in which Un^ni is 
replaced by U {n^) + Urii (n^). Thus, Hdot becomes J2a ^'a^t^^y with the effective levels 
E'^ = Eq + U (n^). Despite of the roughness of HFA, it contains the physics of magnetiza- 
tion due to Coulomb interaction. More important, the approximation allows us to obtain 
a solution including infinite order of tunneling processes, which is crucial for describing the 
Andreev bound states and the supercurrent they carried. However, HFA fails in the Kondo 
regime. As discussed in |T^, if A <C T/< = \/UTe^'^^^°^^^^^ , the spin of the impurity is com- 
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pletely screened by Kondo effect, and there will be no 7r-junction behavior. In this work, 
we will constrain ourselves to the weak coupling case where magnetization effect dominates, 
corresponding to A regime. 

Then most of the formula in the previous section can be transplanted, except the averaged 
occupation number {n„) needs a self-consistent calculation. Notice that (ria) can be derived 
from the retarded Green function of QD as 



i-{ni) = / ^/M[-2ImG^,,(a;)] . 



(18) 
(19) 



Similar to the current formula, (n^) can be divided into two parts, the contribution from 
discrete spectrum and the contribution from continuous spectrum. 
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Since g*" contains unknown quantity (ria) through E'^, the above equations for (no-) should 
be solved self-consistently. 



B. numerical results and discussions 



Next, we present the numerical results for the interacting QD system at zero temperature. 
Fig. 5 shows / vs £^0 and corresponding (n^-) vs £^0 curves for two typical cases, the QD 
without interaction {U = in Fig. 5a) and with strong interaction (f/ ^ F in Fig.Sb). The 
S-QD-S system behaves like one of the above two cases, Depending on the magnitudes of 
the interacting constant U and the coupling strength F, the S-QD-S system may behave 
differently. If C/ < F, the occupation numbers (n-f) and (rij^) are almost equal, leading to the 
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supercurrent always positive with a maximum at Eq = ~U/2, and the system behaves as 
same as the non-interacting one. On the contrast, if f/ > F, a sjTumetry breaking solution of 
Eq.(19) and Eq.(20) is energy preferred, in which (n^) and (nj) are unequal, and QD becomes 
a magnetic dot. Consequently, the effective level occupy a series of configurations 

in the negative plain of Fig.l, so the Josephson current has a small but negative valley in / 
vs Eq curve. We determine the transition border of the above two cases by the parametric 
diagram of / vs {—Eo/r,r/U) in Fig.6. By virtue of electron-hole symmetry, the diagram 
is symmetric to —Eq/U = 0.5. The black area in Fig.6b indicates the range of parameters 
where the S-QD-S system behaves as a vr-junction. One can see from the diagram that the 
system is likely to transfer to vr-junction around Eq = —U /2. 

The temperature dependence are studied in Fig. 7, which shows / vs Eq curves with 
U = 1 for different temperatures. Due to the electron-hole symmetry, only half of the plot 
is shown. The sharp structure at zero temperature is smoothed at finite temperatures, and 
the negative part of the supercurrent vanishes above a critical temperature. These features 
can be understood by taking account of the temperature effect on the supercurrent through 
QD with two spin levels in Fig. 4 and the temperature effect on the averaged QD occupation 
numbers. We also study the diagram of / vs {U,T) where Eq is set to —U/2 (not shown 
here). The transition line from / < to / > can be fitted as U = 0.17 + 7.5T, which is 
consistent with the result derived by non-crossing approximation in |T7 . 



IV. QD IN NON-EQUILIBRIUM DISTRIBUTION 



Now we turn to discuss vr-junction transition caused by non-equilibrium distribution in 
QD. In a recent work by Sun et al. [Q, a mesoscopic four-terminal Josephson junction 
(S-QD-S with two normal leads connected to QD) was studied. By using non-equilibrium 
Green function method, they found that the supercurrent between the two superconducting 
electrodes can be suppressed and even reversed by changing the dc voltage applied across 
the two normal terminals. Here we only take the essential point of that work but omit its 
tedious calculation, by simply assuming that QD has a two-step distribution function as 



i^M = ^[/(^-K) + /(^ + K)]^° 



1 iO < -Vc 

0.5 -Vc<uj < Vc 

UJ>Vr 



(22) 
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corresponding to the limit of r2 = r4 ^ in p4[. Different from p4|, here we allow the two 
spin levels of QD have a Zeeman splitting, i.e., = —Ei = h. The curves of the Josephson 
current J(0 = |) vs the control voltage Vc for different Zeeman splitting h are shown in 
Fig. 8. These curves are step-like because F{lj) is step-like at T = 0, and j{uj) has 5-function 
type discrete spectrum in the range of —A < u < A. Either finite temperature or small 
broadening of Andreev bound states will smooth the curves. 

For h = 0, I reverses its sign around = F, and the magnitude of positive current 
is much larger than that of negative one, which is agree qualitatively with the experiment 



and previous work |]2J]. For h ^ 0, curves of / vs have a peak around Vc = h, with 
the width about 2F, the height about half of that for h = 0. On each side of the peak, 
there is a negative current valley, where the system behaves as vr-j unction. These results 



are also consistent with the calculations for non-equilibrium SFS junction ||2T| , |22[| . One can 
understand these curves by considering the two-step distribution F{u) and CCDOS in Fig.4. 
For example, for the curve of h = 0.2 (marked with "1" in Fig.8), both positive and negative 
Andreev bound states are above the Fermi surface [see CCDOS of {E^,E^) = (0.2, —0.2) ], 
located near h — T and h + T, respectively. When Vc = h, the positive bound state has a 
weight of 0.5, while the negative one has 0, reducing the peak height to one half of that for 
h = 0. For both Vc < h — T ot Vc > h + T, the two Andreev bound states have the same 
weights (either or 0.5), and have little net contribution to the supercurrent. Therefore, 
the negative continuous spectrum of < —A is dominant in the supercurrent, leading to 
the TT-junction transition on both sides of the positive peak. 



V. CONCLUSIONS 



In this work, we have investigated different mechanisms for the vr-junction transition in 
S-QD-S system. From the current formula / = 2sin0 J ^F(uj)j{uj), one can see that the 
vr-junction transition may originate from the change of CCDOS j{uj), or the change of the 
distribution function F{uj), or from the changes of both. The two mechanisms discussed in 
sections II and III, the Zeeman splitting and intra-dot interaction, are involved the change 
of CCDOS j{uj) only. These two mechanisms are closely connected since the intra-dot 
interaction may induce magnetization in QD if the interaction is strong enough. The third 
mechanism studied in section IV involves the change of the distribution function in QD, 
and the interplay of the magnetization with the non-equilibrium distribution in QD. It is 
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interesting that the change of CCDOS and the change of distribution have the similar effect 
on the TT-junction transition, where the positive and negative Andreev bound states cancels 
each other, leaving the negative continuous spectrum dominant in the supercurrent. The 
interplay of the two mechanisms lead to the novel effect that the supercurrent suppressed by 
magnetization of the QD can be partially recovered by a proper non-equilibrium distribution 
of electrons in the QD. 
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VI. APPENDIX 

In this appendix, we discuss the equation of Andreev bound states, i.e. A{u!) — with 
\uj\< A: 

■^c^) - - ^ vfc) h'^' vfe) - S h"- 

Let r = 7A , = EiA , El = — , u = A sin 6' , 9 e (— f , f ), the equation becomes a 
dimensionless form, 

6(e) = (sin^ + 7tane-£i)(sin^ + 7tan^-£2)cos^^ = 7^cos^^ . ((A2)) 

Note that the function oi y — sin^ + 7tane projects Q e (~f > f ) monotonously to 
y e (—00, +00). One can find ©land ©2 in f )> satisfying sinOi + 7tan©i = £1, and 
sin ©2 + 7tan©2 = £2- Suppose ©1 < ©2, we have h{d) ^ for ^ e (-f,©i] U [©2,f), 
and h{Q) < for ^ e (0i,02)- Because 6 (±f ) = 7^, Eq.(A2) has at least two roots 
d\ g(— |,0i]and 02 e [02, f) for 7^ 0. It is straightforward to find the two roots by 
dichotomy method. 

Next, we prove V{B) ^ for < \)[B) < 7^ and 6 e [©2, f ), while b'{e) ^ for 
< b{9) < 7^ and 9 G (— |,0i], so that ^land ^2 are the only two roots in (— |,0i] and 
[02, f ). For 9 e [02, f ), sine + 7tane - £1 ^ 0, and sin^ + 7tane -£2^0. Define 
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X = [(sin^ + 7tan^ — £i)(sin^ + 7tan^ — £2)] , because < h{9) < 7^, so one obviously 
has < X < 7sec^. 

= 2(sine + 7tan e - ^i^)(7 + cos^ 

-2 sin 61 cos d (sin 61 + 7 tan 6' - £1 ) (sin 61 + 7 tan 61 - £2) ((A3)) 
^ 2a;7 — 2 cos 6 ■ = —2 cos 9 ■ x{x — 7 sec 9) ^ . 

Therefore, Eq.(A2) has and only has two roots in the range of ^ e (— |, |) for 7^ 0. 

For = 0, because b' (f ) = -27[2-(£i + £2)], b' (-|) = 27[2 + (£i + £2)], and considering 
the properties of b'{9), one can clearly see that Eq.(A2) has two roots in (— f, f )if 1^1+^21 < 2, 
but only has one root in (— |, |)if |£i + £2] ^ 2. However, this case is irrelevant to the 
Josephson effect. 
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FIGURE CAPTIONS 



Fig. 1 Solution of the equation for Andreev bound states. Parameters are: A = 1, F — 
0.1, = f . (a) and (b) show the two roots Ei and E2 {Ei > E2) vs the configurations 
of QD levels {E^,Ei). (c) and (d) are the diagonal cuts of (a) and (b), showing Ei 
and E2 vs Eq with E^ — E^ = Eq and E^ — —E^ = Eq, respectively. Ei and E2 
can be viewed as two hybrid levels of the electron level oi E — E-^ and the hole level 
of £■ = —El. {El, E2, E^ and E^ are marked as El, E2, Eu and Ed in the plot, 
respectively.) 

Fig. 2 Surface of the Josephson current / vs the configuration of QD levels {E^, E^). Param- 
eters are: A = 1, F = 0.1, = f , T — 0. (a) is the surface graph, while (b), (c) are 
the diagonal cuts, (b) and (c) show / vs Eq with E^ = E^ = Eq and E^ — —E^ = Eq, 
respectively. 

Fig. 3 CCDOS j{u;) for different {E^^Ei) configurations. The row from up to down cor- 
responds to E^ — —0.2, 0, 0.2, and the column from left to right corresponds to 
E^ — —0.2, 0, 0.2. Each of CCDOS contains two types of spectrum: the discrete 
spectrum in the range of |a;| < 1 and the continuous spectrum in the range of |a;| > 1. 
To illustrate discrete spectrum, we broaden 6 functions by 0.01 in the plots. 

Fig. 4 Surfaces of Josephson current / vs (E-^^Ei) configurations at different temperatures: 
T — 10~^, 0.25, 0.50 for (a), (b) and (c), respectively. Other parameters are the same 
as Fig.l. 

Fig. 5 The Josephson current / and the averaged QD occupation number {71^) vs the bare 
QD level £'0, for (a) QD without interaction and (b) QD with strong interaction. 
Parameters are: A = 1, F = 0.1, = f , T = 0, C/ = for (a) and C/ = 1 for (b). 

Fig. 6 Parametric diagram of the Josephson current / vs parameters x = —Eq/U and y = 
T/U. We set A = 1, F = 0.1, = f , T = 0, and change Eq and U to obtain the 
surface graph in (a). The black area of (b) indicates the range of parameter where the 
S-QD-S system behaves as a vr-junction. 

Fig. 7 Temperature dependence of the / vs Eq. Parameters are: A = 1, F = 0.1, (p = 
|, U =1-T = 0.001, 0.05, 0.10, 0.15, and 0.20 for curves marked by 1, 2, 3, 4, and 
5, respectively. 
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Fig. 8 The Josephson current / vs the control voltage Vc for different Zeeman splitting of 
= -El = /i in QD. Parameters are: A = 1, T = 0.1, = f , T = 0; /i = 
0, 0.2, 0.4, 0.6, and 0.8 for curves marked with 0, 1, 2, 3, and 4, respectively. 
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